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Abstract. We report the results of a high temporal resolution spectroscopic monitoring of the flare star AD Leo. 
During 4 nights, more than 600 spectra were taken in the optical range using the Isaac Newton Telescope (INT) 
and the Intermediate Dispersion Spectrograph (IDS). We have observed a large number of short and weak flares 
occurring very frequently (flare activity > 0.71 hours"^). This is in favour of the very important role that flares 
can play in stellar coronal heating. The detected flares are non white-light flares and, though most of solar flares 
belong to this kind, very few such events had been previously observed on stars. The behaviour of different 
chromospheric lines (B aimer series from Hq to Hn, Call H & K, Nai Di & D2, Hei 4026 A and Hei D3) has 
been studied in detail for a total of 14 flares. We have also estimated the physical parameters of the flaring plasma 
by using a procedure which assumes a simplifled slab model of flares. All the obtained physical parameters are 
consistent with previously derived values for stellar flares, and the areas - less than 2.3 % of the stellar surface - 
are comparable with the size inferred for other solar and stellar flares. Finally, we have studied the relationships 
between the physical parameters and the area, duration, maximum flux and energy released during the detected 
flares. 
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1. Introduction 

Stellar flares are events where a large amount of en- 
ergy is released in a short interval of time, taking place 
changes at almost all frequencies in the electromagnetic 
spectrum. Flares are believed to be the result of the re- 
lease of part of the magnetic energy stored in t he corona 
through magnetic reconnection (see reviews by iMirzovanI 
Il984t iHaisch et ai]ll99lt lGarcia-Alvare3l200o() . However, 
the exact mechanisms leading to the energy release and 
subsequent excitation of various emission features remain 
poorly understood. Many types of cool stars produce 
flares, sometimes at levels several orders of magnitude 
more energetic than their solar counterparts l)PetterserJ 
ll989HGarcia- Alvarez et"al]l2002l) . In dMe stars (UV Ceti- 
type stars) optical flares are a common phenomenon. On 
the contrary, in more luminous stars flares are usually only 



detec ted through UV or X-ray observations l|Dovle et alJ 
though some optical flares have been observed 



in young early K 


dwarfs like LQ Hya 


and PW And 


l|Montes et al., 


1999,: 


,L6dcz- Santiago et al. 200^). 
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One would like to trace all the energetic processes 
in a flare to a common origin, though the released en- 
ergy can be very different . The largest so lar flares in- 
volve energies of 10 ergs (lGershbergili989[l . Large flares 
on dMe stars can be t wo or d ers of magnitude large r 
l|Dov1e fc Mathiouda,kiFl ES^S iBvrne McKa^ ^2^, 
while very energetic flares are produced by RS CVn bi- 
nary systems, where the total energy may exceed 10^^ ergs 
llDovle et al.ll992tlrbhig et al.ll99llGarcia-Alvarez et alJ 
1200.'^ . Such a change in the star's radiation field modi- 
fies drastically the atmospheric properties over large areas, 
from p hotospheric to coronal layers. Models bv lHoudebinj 
11X99211 indicate that heating may be propagated down 
to low photospheric levels, with densities higher than 
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Table 1. Observing log INT/IDS (2-5 April 2001). 



R1200B R1200Y 

Night N UT toxp(s) SNR N UT toxp(s) SNR 

start-end min— max min— max start-end min— max min-max 















H/3 


H7 


H(5 


Can H 
Can K 


Hs 

Hg 

Hit, 








Hq 


Nal Di 
Nal D2 
He I D3 


1 


18 


00:02- 


-01:47 


180 


-300 


76-85 


53-58 


45-49 


33-37 


28-31 













2 


91 


20:19- 


-01:58 


60- 


120 


46-73 


32-48 


26-40 


20-29 


17-25 













3 


246 


20:10- 


-02:37 


15- 


120 


33-89 


21-58 


18-51 


13-38 


12-31 













4 


104 


20:28- 


-22:44 


15- 


120 


27-72 


18-48 


15-41 


11-30 


10-25 


189 


23:38-03:13 


5-120 


32-172 


23-114 



10 cm However, electrons with energies in the MeV 
range would be required to attain such depth. 

Spectral emission lines are the most appropriate di- 
agnostics to constrain the physical propertie s and mo- 
tions of flaring plasmas (see iHoudebind EoOSl and refer- 
ences therein). Unfortunately, only a few sets of observa- 
tions with adequate time and spectral resolution are so far 

sL 

un- ^ „ 

2002t iHawlev et~ai.l l200,1^ . Furthermore, there are also 
few attempts t o deriv e the physical parameters of the 
flare plasmas llDonati -Falchi et al. 1985; Kunkel 197(| 
CTershberd ll 974 lKa,tsoval ll 99(11: I.Tevremovic et a,]lll998ll . 
Garcia- Alvarez et al.l l|2002|) obtained the first detailed 



tions with adequate time and spectral resolution are so tar 
available teodo no et 'aI]ll989HHawlev fc Pette rsenlll99lt 
Houdcbinc 19 92,:lGiiiin et alJll994bl: lGarcfa- Al varez et alJ 



trace of physical par ameters during a l argo optical flare, 
on AT Mic, using the I Jevremovic et al.l l)l998l) procedure. 
Good quality spectroscopic observations, as those ana- 
lyzed in this work for AD Leo, and a correct description 
of the different flare components are required to constrain 
numerical simulations. 

AD Leo is well-known for being a frequent source 
of flares. It has been subject of numerous studies in 
the optical, EUV and X-rays because of its nature 
as on e of th e most active M d w arfs llPettersen et al, 
1984 Il990t iLang fc Willsor] |l98(it iGiidel et al 
1989 : iHawlev fc Pettersen] " 
20031 lAbada-Simon et al " 



11991b 
I1997 



Hawlej^^^^ 



iGullv et alJ 1997; 



Favata et alJ I2OO0I: ISanz-Forcada fc Micelaf 
van den Besselaar et all l2003t iMagg o et alJ 



1995 



2002; 



12004 



va.n rien hes.se I aar eta 1. 1 l^( Ml.it IMagg 

Bobrade fc SchmittI Mi^ ISmith et a,lJ 12005^. AD Leo 
(GJ 388) is classified as dM3Ve llHenrv et al.lll994) . This 



star has an unseen companion - detected using speckle 
interferometry - with a period of a bout 27 years which 
is expected to have a very low mass i Balega et alJ ll9841. 
AD Leo is located in the immediate solar neighbourhood, 
at a distance of ^ 4.9 pc (from the ground-based parallax 
ijGrav fc Johansonlll99l)) - it was not a Hipparcos target). 
Its high activity is probably due to its high ro tation rate 
(-Pphot ~ 2.7 days, ISpiesman fc HawlevI I198(tI) . In fact, 
its rotational velocity (?;sim = 6.2 ± 0.8 km s~^) places 
AD Leo among the tail of rare fast rotating M dwarfs 
jDelfosse et al.lll998l) . The radius and mass of A D Leo are 
R ~ O.44i?0 and M ~ 0.40Mq, respectively (|Pettersenl 



ll976tlF^vata et ai]l2000 l). At this mass, stars are expected 
to have a substantial radiativ e core, so that the interio r 
structure is st ill "sola r-like" l|Ghabrier fc Baraffel 
ISaa,r fc Linsk^ (1985) detected, through infrared line 
measurements, photospheric fields on AD Leo showing 
the presence of strong magnetic fields. They inferred that 
73 % of the AD Leo's surface is covered by active regions 
outside of dark spots containing a mean field strength of 
B = 3800 ± 260 G. 

Here we present the results of a long spectroscopic 
monitoring of AD Leo that was carried out using an inter- 
mediate dispersion spectrograph and high temporal reso- 
lution. This work has considerably extended the existing 
sample of stellar flares analyzed with good quality spec- 
troscopy in the optical range. Details about the techni- 
cal information of the observations and data reduction 
are given in § |2| Section O describes the analysis of the 
observations and the detected flares, including equivalent 
widths, line fluxes, released energy, line profiles and asym- 
metries. In § 0] we present the main physical plasma pa- 
rameters obt ained for the obse r ved fl ares using the code 
developed bv I Jevremovic et al.l l)l998l) . The discussion of 
the results and conclusions are given in § El where the 
relationships between the physical and observational pa- 
rameters are also analyzed. The preliminary result s for 
this st ar a nd V1054 Oph were pr esented bv .Montes et alJ 
l|2003() and lCrespo-Chacon et all l(2004() . 



2. Observations and data reduction 

The data were taken during the MUlti-SIte Continuous 
Spectroscopy (MUSICOS) 2001 campaign. It involved 
observations at two sites: El Roque de los Muchachos 
Observatory from La Palma (Spain) and SAAO (South 
Africa). However, due to poor weather conditions, the 
only useful data were those taken at La Palma. This 
observing run was carried out with the 2.5 m Isaac 
Newton Telescope (INT) during 2-5 April 2001. The 
Intermediate Dispersion Spectrograph (IDS) was utilized 
together with the 2148x4200 EEVlOa CCD detector. Two 
gratings were used: R1200B (every night) and R1200Y 
(second half of the last night). The wavelength covered by 
R1200B (blue spectrum) ranges from 3554 A to 5176 A 
(including the Balmer lines from H/3 to Hn as well as 
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Reference star: Gl 687B (M3.5V) ^ 
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Fig. 1. Observed spectrum of AD Leo at the maximum of the strongest flare detected with the R1200B grating (flare 
2, see S 13.1(1 and in its quiescent state. The observed spectrum of the reference star Gl 687B is shown at the bottom. 
The chromospheric lines are identified. 
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Fig. 2. As Fig. ^but for the spectra taken with the R1200Y grating, the strongest flare detected using this spectral 
configuration (flare 12, see S 13. 1(1 . and the reference star GJ 725B. 



the Can H & K and Hei 4026 A lines). The recip- 
rocal dispersion of these spectra is 0.48 A/pixel. The 
wavelength covered by R1200Y (red spectrum) ranges 
from 5527 A to 7137 A (including the Ha, Nai Di & 
D2 and He I D3 lines), being the reciprocal dispersion 
0.47 A/pixel. The spectral resolution, determined as the 
full width at half maximum (FWHM) of the arc compari- 



son lines, is 1.22 A for the blue region and 1.13 A for the 
red one. 

We took series of spectra with short exposure times: 
from 15 to 300 s for R1200B and from 5 to 120 s for 
R1200Y. The spectra in each series were separated only 
by the CCD readout time (less than 60 s) in order to ob- 
tain the highest temporal resolution as possible. During 
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Fig. 3. Temporal evolution of the EW of the H/3 line. The strongest observed flares are labeled. 
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Fig. 4. As Fig.Obut for the Ha line. The observed flares 
have been numbered following the order in Fig. |2| 



4 nights, a total of 459 spectra of AD Leo were ob- 
tained with R1200B and 189 with R1200Y. The observ- 
ing log (Table ^ lists: the number of spectra taken ev- 
ery night [N], the universal time (UT) at the beginning 
and end of the AD Leo observations, the minimum and 
maximum exposure time (texp) of the spectra, and the 
signal-to-noise ratio {SNR) of the continuum near each 
line region. These quantities are given for both gratings 
(R1200B and R1200Y). Note that, despite the short expo- 
sure times, the SNR is large enough to perform a reliable 
analysis. 



The reduction has been done following the standard 
procedure: bias and dark subtraction, flat-field correction 
using exposures of a tungsten lamp, cosmic rays correc- 
tion, and sky background subtraction from the region of 
the aperture chosen for doing an optimal extraction of 
each spectrum. The software packages of IRAF^ have been 
used. The wavelength calibration has been done by using 
spectra of Cu-Ar lamps. Finally, all the spectra have been 
normalized to their maximum flux value of the observed 
continuum. 

3. Analysis of the observations 

In Fig. ^ we have plotted the observed blue spectrum of 
AD Leo in its quiescent state (minimum observed emis- 
sion level) and at the maximum of the strongest flare de- 
tected with the R1200B grating (flare 2, see §12111. Fig. El 
is similar to Fig. ^ but for the red spectra, taken with 
the R1200Y grating, and the strongest flare detected with 
this spectral configuration (flare 12, see S I3.1|I . In addition, 
the spectra of Gl 687B and GJ 725B have been plotted as 
examples of inactive stars with a spectral type and lu- 
minosity class (M3.5V) very similar to those of AD Leo 
(M3V). Note that the origin of the Y-axis is different for 
every spectrum in order to avoid the overlapping between 
them. 

AD Leo shows a high level of chromospheric activity. A 
strong emission in the Balmer series and the Can H & K 
lines is observed even in its quiescent state (see Fig. ^ 
and El . The greatest emission in these lines is found at 
the maximum of the detected flares. The Hei D3 and 

^ TRAP is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Table 2. Minimum and maximum relative error in the EW of the different chromospheric lines. 



5 



Night A_BW,oi (%) 

(min-max) 



H/3 H7 Hi5 HciA4026A CaiiK CaiiH+Ht Hg Hg Hio Hn Ha NaiDi Nai D2 HCID3 



1 10-13 6-9 8-14 20-46 5-8 5-8 9-20 9-17 11-24 16-37 - - - - 

2 9-15 6-10 7-17 19-57 5-9 5-10 9-23 9-22 12-34 17-62 - - - - 

3 10-17 6-12 8-20 24-109 5-13 5-13 10-28 9-29 12-41 17-61 - - - - 

4 10-20 6-15 8-23 27-152 5-14 6-15 10-34 10-39 12-56 18-145 9-20 2-5 1-5 32-328 



He I 4026 A lines also show emission above the contin- 
uum, which is more noticeable during the strongest flares. 
The Nai Di & D2 absorption lines only present a slight 
filled-in profile in the quiescent state, though a low emis- 
sion is observed in their core during flares. High resolu- 
tion spectroscopic observations of the quiescent state of 
AD Leo confirm the presence of chromospheric emission 
lines that range from a weak emission in the center of 
the Can IRT and the Nai Di & D2 lines to a strong 
emission in the Hei D.':;, Mgll h fc k, Call H fc K and 
H I Balmcr l ines JPettersen fc Cole man"l98l': 'Dovle'igST"; 
ICresp o- C hacon et al.ll2006|) . In addition, Sundland et al. 
I|1988j) suggested that Balmer lines are the most impor- 
tant components to chromospheric radiation loss. 

AD L eo is well-known for having strong flares (see, for 
example. iHawlev fc PettersenlllflQlj) . However, after over- 
lapping the normalized spectra and comparing the depth 
of the absorption lines and the shape of the continuum, 
we have not detected any noticeable continuum change 
(see also Fig. and EJ. Therefore the flares analyzed in 
this work are non white-light flares and, within the stud- 
ied wavelength range, they only affect the emission in the 
chromospheric lines. Despite this kind of flares is the most 
typical in the Sun, in which no detectable signature in 
white-light is generally observed, very few such events 
had been previously detected on stars (see iButler et al.l 
Ll986: Houdebine 2Q03i) . Conversely to the Sun, stellar non 
white-light flares have been seldom observed because very 
little time has been dedicated to spectroscopy in compari- 
son to photometry. However. .Houdebine 11992) suggested 
that the low frequency observed for stellar non white-light 
flares could be due to a strong contrast effect. In other 
words, many flares have been detected as white-light flares 
on dMe stars because of the relative weakness of their pho- 
tospheric background, but they would have been classified 
as non white-light flares on the Sun because of its much 
higher photospheric background. 

3.1. Equivalent widths 

The equivalent width {EW) of the observed chromo- 
spheric lines has been measured very carefully in order 
to detect possible weak flares in our observations. Two 
methods have been tried out. The first one consists of 
obtaining the EW with the routine SPLOT included in 
IRAF, taking exactly the same wavelength limits for each 



line in all the spectra. The second method uses the routine 
SBANDS of IRAF, which generally introduces less noise 
in the EW measurement because of taking into account a 
larger number of points to calculate the continuum value. 
For this reason, we only present the results obtained with 
SBANDS. Note that the line regions have been taken wide 
enough to include the enhancement of the line wings dur- 
ing flares. 

SBANDS uses Eq.Ulto obtain the EW of a hue; where 
W\ is the width at the base of the line, d is the reciprocal 
dispersion (pixel size in A), i^c.i is the flux per pixel in the 
continuum under the line, is the observed flux in the 
pixel i of the line, and n is the number of pixels within the 
line region. 



i^l 



(1) 



The EW uncertainty has been estimated using Eq. 13 
We have considered AWi = A and Ad = 0.01 A/pixel. 
Eq. 121 has been obtained by applying the stan- 
dard quadratic error propagation theory to Eq. ^ 
taking into account that n = W\/d and assuming 
AFi,i « AFe,i = Fe,i/SNR. 



(2) 



Ci = 
C2 = 
C3 = 

C4 = 



where 
dEW 



dWi 



dEW . , EW-Wi . . 

dEW . ^ Wi- EW 
A i<c,i = 



dEW 
dFiA 



A F14 = 



SNR 



SNR 



Table |21 shows the minimum and maximum values of 
the relative error in EW {AEWrd) for each chromospheric 
line and every night. The AEW,-e\ of the Can H fc K, 
Nai Di fc D2, H/3, H7 and US lines is between 1 % and 
20 %. For Hg and Hg the AEW,c\ is less than 40 %, while 
for Hio and Hn is higher than the uncertainty estimated 
for the other Balmer lines (up to 50 % in the case of Hio 
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Table 3. JD of the flare onset {JDsta-ct): total duration of the detected flares and length of their phases (seen in the 
H/3 line in the case of flares 1 to 11 and seen in the Ha line in the case of flares 12 to 14). The time delay between 
the maximum of the chromospheric lines and the one of H/3 or Ha, depending on the spectral configuration, is also 
shown. Blanks are given when the beginning, maximum and/or end of the flare was not observed in the line under 
consideration. 



Flare 


J-Dstait (days) 




Duration (min) 


Delay of the 


maximum (min) 




(2452000+) 














Total 


Impulsive phase 


Gradual decay 


Can K 


Hoi 4026 A 






(H/3) 


(H/3) 


(H/3) 






1 


2.513±0.003 




13±7 






0+6 


2 


3.3903±0.0016 


25±4 


11±5 


14+4 


0+5 


0+5 


3 


3.5321±0.0016 




10±4 




0+4 


0+4 


4 


3.5608±0.0013 


31±3 


10±4 


22+3 


0+4 


0+4 


5 


4.3781±0.0010 


18±4 


6±2 


11+3 


0+1 




6 


4.4193±0.0005 


22±1 


4±2 


17+2 


5+3 


0+3 


7 


4.4804±0.0010 


25±2 


9±3 


16+2 


3+2 


1+2 


8 


4.5479±0.0009 




10±2 




1+1 


1+1 


9 


4.5978±0.0005 


17±2 


5±1 


12+2 


2+1 




10 
11 


5.3968±0.0005 


14±1 


1±1 


12+1 


3+1 




Flare 


J-Dstart (days) 




Duration (min) 


Delay of the 


maximum (min) 




(2452000+) 
















Total 


Impulsive phase 


Gradual decay 


Nai Di, D2 Hoi D3 






(Ha) 


(Ha) 


(Ha) 






12 


5.5341±0.0003 


31±2 


5.9+0.9 


25+2 


0.9+0.9 




13 


5.5763±0.0005 


14±1 


2+1 


12+1 


1 + 1 


1 + 1 


14 















and even larger than 100 % in the case of Hn). In addi- 
tion, the /\EW^c\ of He i D3 and He i 4026 A is frequently 
greater than 100 %. For this reason, the results obtained 
for the He i lines are less reliable. Besides, because of the 
short exposure times and the relative faintness of AD Leo 
in the blue, the region blueward Hio is usually too noisy 
for reliable measurements. 

The observed flares have been detected using the H/3 
and Ha lines (for the blue and red spectra, respectively). 
H/3 and Ha have been chosen because they suffer from 
large variations during flares and are the chromospheric 
lines with the best SNR in each spectral configuration. 

Fig.|2|shows the temporal evolution found for the EW 
of the H/3 line. The observed fiares are marked with num- 
bers. Other smaller changes can also be seen. The tem- 
poral evolution found for the EW of the Ha line is given 
in Fig. 01 No strong variations have been observed during 
our spectroscopic monitoring. However, 14 short and weak 
flares have been detected: 11 among the blue spectra and 
3 among the red ones. The Julian Date [JD) at the be- 
ginning of each flare ( JUstart) is shown in Table |21 These 
flares have been also observed in the other chromospheric 
lines. Nevertheless, it is more difficult to distinguish be- 
tween flares and noise for lower SNR and/or weaker lines. 
In addition, a modulation with a period of ~ 2 days is also 
observed in the quiescent emission of AD Leo, though the 
photometric period of this star was found to be 2.7 days 
I Spiesman fc Hawlevl[l98^ . 

We want to emphasize that the flare frequency is very 
high, and sometimes new flares take place when others are 
still active (see, for example, flares 8 and 11 and the points 



before flares 2 and 7 in Fig. El . It is also remarkable that 
the flares in night 4 seem to erupt over the gradual decay 
phase of other stronger flare (Fig.^. For calculating the 
flare frequency we have taken the total observed time as 
the sum of UTend - UTgtart for all the nights (see Table^) . 
Note that, within a night, the series of observations are not 
continuous in time. Therefore the flare frequency should 
be greater than the obtained one. This implies that the 
flare activity of AD Leo, found from the variations of 
the chromospheric emission lines, is > 0.71 flares/hour. 
This value is slightly higher than those that other au- 
thors measured by using photometry: Moffe^(^3) ob- 
served 0.42 flares/hour, IPettersen et al.l lll984tl detected 
0.57 flares/hour, and iKonstantinova-Antova fc Antovl 
(1995) found 0.33 - 0.70 flares/hour. Given that the de- 
tected flares are non white-light flares, we conclude that 
non white-light flares may be more frequent on stars than 
white-light flares, as observed in the Sun. In fact, as the 
total time of the gaps between the series of observations is 
similar to the total real observed time, the flare activity of 
AD Leo could even be the double of that obtained. As far 
as we know, this is the first time that such a high fiare fre- 
quency is inferred from the variation of the chromospheric 
emission lines. 

3.2. Equivalent widths relative to the quiescent state 

In order to compare the behaviour of the different chromo- 
spheric lines, we have used the equivalent width relative 
to the quiescent state {EWRQ)^ defined as the ratio of 
the EW to the EW in the quiescent state. Fig. [S] and 
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Table 4. EWRQ at flare maximum (-EVF-RQmax) for the different chromospheric lines. Blanks are given when the 
flare maximum was not detected in the line under consideration. Flares with no maximum data in any line have been 
omitted. 



Flare EWRQ, 





H/3 


H7 


H<5 


Hs 


H9 


Hio 


Hii 


Can H + He 


Can K 


He I 4026 A 


1 


1.26±0.18 


1.28±0.11 


1.47±0.17 


1.53±0.22 


1.52±0.20 


1.7±0.3 


1.6±0.4 


1.16±0.09 




1.5±0.5 


2 


1.65±0.25 


1.52±0.14 


2.0±0.3 


2.0±0.3 


1.8±0.3 


1.5±0.3 


2.0±0.7 


1.30±0.11 


1.14±0.09 


1.7±0.6 


3 


1.09±0.19 


1.06±0.10 


1.10±0.17 


1.11±0.21 


1.13±0.22 


0.99±0.23 


1.0±0.4 


1.1±0.1 


1.01±0.09 


1.3±0.5 


4 


1.23±0.20 


1.25±0.12 


1.44±0.21 


1.5±0.3 


1.5±0.3 


1.6±0.3 


1.6±0.6 


1.13±0.10 


1.06±0.09 


1.2±0.5 


5 


1.08±0.20 


1.08±0.13 


1.14±0.21 


1.1±0.3 


1.1±0.3 


1.0±0.3 


1.5±0.6 


1.06± 0.13 


1.01±0.12 




6 


1.4±0.3 


1.34±0.17 


1.6±0.3 


1.8±0.5 


1.4±0.4 


1.5±0.5 


1.3±0.5 


1.23±0.16 


1.14±0.13 


1.6±0.9 


7 


1.32±0.24 


1.32±0.17 


1.5±0.3 


1.4±0.4 


1.7±0.4 


1.5±0.5 


1.9±0.8 


1.15±0.15 


1.09±0.14 


1.9±1.1 


8 


1.27±0.23 


1.25±0.14 


1.49±0.25 


1.5±0.3 


1.5±0.3 


1.4±0.4 


1.5±0.6 


1.16±0.13 


1.06±0.11 


1.6±0.8 


9 


1.16±0.21 


1.07±0.13 


1.18±0.21 


1.2±0.3 


1.3±0.3 


1.1±0.3 


1.0±0.4 


1.13±0.14 


1.06±0.14 




11 


1.12±0.22 


l.lOitO.15 


1.19±0.24 


1.3±0.3 


1.1±0.3 


1.2±0.4 


1.0±0.5 


1.05±0.13 


1.07±0.13 




Flare 




EWRQ 




















Ha 


Nal Di 


Nai D2 


He I D3 














12 


1.09±0.15 


0.95±0.04 


0.95±0.03 
















13 


1.13±0.18 


0.98±0.04 


0.98±0.04 


1.5±1.1 















show the EWRQ of several lines for a representative sam- 
ple of the detected flares (those in which all the phases 
have been observed). The chromospheric lines plotted in 
these two figures are those with less uncertainty in the 
EW measurement. We can observe two different kinds of 
flares: for some of them (sec flares 6, 9, 11, 12 and 13) the 
gradual decay phase of the Balmer lines is much longer 
than the impulsive phase; in contrast, other flares (2, 4, 5 
and 7) are less impulsive. This resembles the classiflcation 
of solar flares made bv iPallavicini et al. (1977): eruptive 
flares or long-decay events and confined or compact flares. 
However, all these flares always show the same behaviour 
in the Can H & K lines, that is a slow evolution through- 
out all the event. The evolution of the Nai Di & D2, 
He I D3 and Hei 4026 A lines is quite similar to that of 
the Balmer series. However, it is less clear for the Hei lines 
because of having a very high error sometimes. 

In some flares (5, 6, 11, 12 and 13) we can see two 
maxima in the Nai and Balmer lines, but it is not so ev- 
ident for the Hei and Can lines. Sometimes variations 
at shorter time scales are also detected during the grad- 
ual decay phase, in which different peaks, decreasing in 
intensity, are observed in the EWRQ (see, for instance, 
flare 12). This could be interpreted as the succession of 
different reconnection processes, decreasing in efficiency, 
w ithin the same flare - foll owing; the original suggestion 
of IKodd & Pneuma - or as a series of flares orig- 

inated as a result of wave disturbances coming along the 
stellar surface from the first flare region. 

Table |3I contains: the total duration of the detected 
flares, the length of their impulsive and gradual decay 
phases, and the time delay between the maximum of each 
chromospheric line and that of H/3 or Ha (depending on 
the spectral configuration). The detected flares last from 
14 ± 1 up to 31 ± 3 min. Regarding the first maximum 
of emission, we have found that the Balmer series reach it 
simultaneously while the rest of the lines are delayed. The 
delay is negligible for the Hei and Nai lines 1 min) 



but it is very evident for Can H & K (up to 5 ± 3 min). 
It seems that this delay is greater for the flares with a 
shorter impulsive phase of the Balmer lines. However, the 
total duration of the flare does not seem to be related. For 
the first five flares, we cannot assume that the delay of the 
maximum is completely zero because their uncertainties, 
particularly those of Hei 4026 A, are even greater than 
the delay found in the other cases. The moment in which a 
line reaches its maximum is related to the height where the 
line is formed above the stellar surface. This is also related 
to the temperature that characterizes the formation of the 
line. According to the line formation models in stellar at- 
mospheres, the Can H & K lines are formed at deeper 
and cooler layers than the Balmer series. Therefore, the 
gas that is heated and evaporated into the ne wlv formed 
loop after magnetic reco nnection l)Cargill fc Pries t 19*83 
iForbes fc Malherbdll986^ cools and reaches the formation 
temperature of the Balmer series before the one of the 
Can H & K fines. Houdebin c (2003) found that the rise 
and decay times in the Call K and H7 lines obey good 
relationships, which implies that there is a well-defined un- 
derlying mechanism responsible for the flux time profiles 
in these lines. 



Table 0] lists the EWRQ of the chromospheric lines 
at their maximum in each flare {EW RQ-aiayi)- The lower 
the wavelength, the greater EWRQmax is found for the 
Balmer lines. However, H/3 and H7 have a very similar 
variation. For Hei 4026 A and Hei D3 the EWRQma.^ 
is analogous to that of the Balmer series, whereas for the 
Can H & K and Nai Di & D2 lines is smaller. Our results 
show that the duration of flares tends to be larger when 
the EWRQraax of the Balmer lines is greater. This effect is 
more noticeable for the less impulsive flares. Nevertheless, 
there are no clear relationships between the EWRQ^ax 
of the other lines and the duration of the flare. 
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Fig. 5. EWRQ of the flares (detected using the R1200B grating) in which all the phases have been observed (lines: 
H/3, H7, m, Hg, Can H + He, Can K, Hei 4026 A). 
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Fig. 6. As Fig.|31but for the flares detected using the R1200Y grating (lines: Ha, Nai Di, Nai D2, Hei D3). 



3.3. Line fluxes and released energy 

In order to estimate the flare energy released in the ob- 
served chromospheric lines, we have converted the EW 
into absolute surface fluxes and luminosities. 

The absolute line fluxes (F) have been computed from 
the EW measured for each line and its local contin- 
uum. The absolute flux of the continuum near each line 
ha s been determined making use of the method given 
bv IPettersen fc HawlevI 1)19891) . For AD Leo, they inter- 
polate between the R and I filter photometric fluxes to 
estimate the observed flux in the continuum near 8850 A. 
This value is then transformed into a continuum sur- 
face flux of i^(8850 A) = 6.5 x 10^ erg s-^cm-^A"!, us- 
ing 0.44 Rq for the radius and 0.203 arcsec for the par- 
allax. Direct scaling from the flux-calibra t ed spe ctrum of 
AD Leo, showed bv IPettersen fc HawlevI l)l989|) . has al- 
lowed us to estimate the continuum surface fluxes near 
the emission lines of interest. We multiply this continuum 
value by the measured EW to determine the line surface 
flux. The absolute fluxes of the Balmer lines at the dif- 
ferent flare maxima (^max) are listed in Table |3] These 
values will be used in § 0] to calculate the Balmer decre- 
ments. Note that each F^ax only shows the contribution 
of the flare to the surface flux (the contribution of the qui- 
escent state has been subtracted). Since our aim is to iso- 
late the contribution of the flares from the total emission, 
we have chosen the quiescent state for each of them as a 
pseudo-quiescent, that is the minimum emission level just 
before or after the flare. In this way, variations originated 
by other processes, which can contribute to the observed 
emission, are not taken into account. Even if these pro- 
cesses were magnetic reconnections within the same active 
region, assuming that the solar flare model is also applica- 
ble to the stellar case, the new emi tting plasma would be 
placed in a different post- fiare loop llCargill fc Priestll983t 
iForbes fc Malherbelll984 . Therefore, in first approxima- 
tion, its physical properties (see § could be separately 
studied. 



We have also converted the absolute fluxes into lu- 
minosities (L). Table [Sj lists the contribution of every 
flare to the energy released in the Balmer lines. The re- 
leased energy has been calculated by numerically integrat- 
ing the luminosity from the beginning to the end of the 
flare and subtracting the contribution of its correspond- 
ing pseudo-quiescent state. The energy released during the 
impulsive and gradual decay phases, which has been cal- 
culated in the same way, is given as well. 

3.4. Line profiles and asymmetries 

We focus this section on the Balmer series and Can H 
fc K lines. The intermediate spectral resolution of the ob- 
servations has not allowed us to study the profile of the 
He I and Na l lines because of their low intensity. 

We have plotted the temporal evolution of the ob- 
served line profiles during the detected flares (see the ex- 
amples given in Fig. No evidence for line-shifts has 
been found. However, the reported flares may be too weak 
for measuring line-shifts using the available spectral reso- 
lution (see § O . It seems that flares affect the core of the 
lines before the wings. The emission in both of them (core 
and wings) increases during the impulsive phase, shows 
the largest value when the EWRQ of the line reaches the 
maximum, and decreases during the gradual decay. The 
emission in the wings of the Ca 11 H fc K lines hardly rises 
during the observed flares, being only noticeable during 
the strongest ones (see Fig.[7|). In fact, their broadening is 
negligible even at flare maximum. If the flare is too weak, 
as the flare 5, the profile of the Balmer lines shows the 
same behaviour as that found for the Can H fc K lines. 
In contrast, when the released energy is high enough to 
produce changes in the wings of the Balmer series (flares 
2, 6 and 7), the core of these lines seems to rise up to a 
constant value (which is independent of the EWRQmax) 
whereas the emission in the wings depends on the flare 
energy. Table lists the width at the base of the lines 
of interest during the quiescent state and the maximum 
of the strongest observed flare (flare 2) . The fact that the 
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Table 5. Surface fluxes of the Balmer lines at flare maximum (i^max)- The contribution of the quiescent state has 
been subtracted. Flares with no maximum data in any line have been omitted. 



Flare 








ax (10^ 


org s 








Hq 


H/3 


H7 


H5 


H8 


H9 


Hio 


Hii 


1 




0.90 


0.64 


0.56 


0.29 


0.23 


0.109 


0.070 


2 




1.82 


1.10 


0.90 


0.47 


0.33 


0.108 


0.105 


3 




0.34 


0.19 


0.15 


0.08 


0.053 


0.011 


0.013 


4 




0.70 


0.52 


0.41 


0.22 


0.17 


0.120 


0.079 


5 




0.34 


0.23 


0.19 


0.10 


0.04 


0.079 


0.061 


6 




1.19 


0.70 


0.62 


0.32 


0.16 


0.084 


0.040 


7 




0.94 


0.62 


0.43 


0.19 


0.25 


0.076 


0.090 


8 




0.74 


0.48 


0.42 


0.16 


0.14 


0.072 


0.053 


9 




0.57 


0.27 


0.30 


0.11 


0.11 


0.070 


0.036 


11 




0.35 


0.07 


0.18 


0.08 


0.03 


0.050 


0.010 


12 


1.10 
















13 


1.40 

















Table 6. Energy released in the Balmer lines during every flare phase. The total released energy is also shown. The 
contribution of the quiescent state has been subtracted. We have omitted the flares with no data about their beginning, 
maximum and end. 



Flare Phase Released energy (10^^ erg) 







Hq 


H/3 


H7 


H<5 


He 


H9 


Hio 


Hii 


1 


Impulsive 




2.5 


1.2 


1.6 


0.8 


0.7 


0.3 


0.2 




Gradual 




















Total 


















2 


Impulsive 




3.3 


2.3 


1.4 


0.9 


0.9 


0.2 


0.3 




Gradual 




4.9 


3.1 


2.0 


1.2 


1.2 


0.4 


0.5 




Total 




8.2 


5.4 


3.4 


2.2 


2.1 


0.6 


0.7 


3 


Impulsive 




1.2 


0.7 


0.7 


0.3 


0.2 


0.1 


0.0 




Gradual 




















Total 


















4 


Impulsive 




1.9 


1.3 


0.9 


0.6 


0.5 


0.4 


0.3 




Gradual 




2.2 


1.6 


0.6 


0.9 


0.5 


0.5 


0.5 




Total 




4.1 


3.0 


1.5 


1.5 


1.1 


0.9 


0.8 


5 


Impulsive 




0.6 


0.2 


0.3 


0.2 


0.1 


0.3 


0.2 




Gradual 




1.0 


0.6 


0.4 


0.2 


0.1 


0.4 


0.2 




Total 




1.5 


0.8 


0.7 


0.4 


0.2 


0.7 


0.4 


6 


Impulsive 




1.5 


0.7 


0.9 


0.5 


0.2 


0.1 


0.0 




Gradual 




6.2 


3.3 


3.1 


1.8 


1.0 


0.3 


0.1 




Total 




7.7 


3.9 


4.0 


2.3 


1.2 


0.4 


0.1 


7 


Impulsive 




1.8 


1.2 


1.1 


0.4 


0.5 


0.2 


0.2 




Gradual 




5.4 


3.4 


3.2 


1.3 


0.8 


0.6 


0.6 




Total 




7.2 


4.6 


4.4 


1.7 


1.2 


0.8 


0.9 


8 


Impulsive 




2.0 


1.3 


1.0 


0.4 


0.4 


0.2 


0.1 




Gradual 




















Total 


















9 


Impulsive 




0.6 


0.4 


0.3 


0.0 


0.0 


0.1 


0.0 




Gradual 




1.5 


0.9 


0.7 


0.1 


0.1 


0.3 


0.2 




Total 




2.1 


1.3 


1.0 


0.1 


0.1 


0.5 


0.2 


11 


Impulsive 




0.2 


0.1 


0.1 


0.0 


0.0 


0.0 


0.0 




Gradual 




1.8 


0.6 


1.2 


0.2 


0.2 


0.4 


0.1 




Total 




2.0 


0.7 


1.3 


0.2 


0.2 


0.5 


0.1 


12 


Impulsive 


3.1 


















Gradual 


10.3 


















Total 


13.4 
















13 


Impulsive 


1.0 


















Gradual 


6.5 


















Total 


7.5 

















width of the Ca ii lines is much less sensitive to flares than 
that of the Balmer lines sug gests that the broadening ma y 
be due to the Stark effect l(Bvrnelll98cH iRobinsonlllQSfll) . 

Nevertheless, mass motions can also be present (see be- 
low). 

Fig.lHlshows, as an example, the bisector of the H/3 and 
Can K lines during the flare 2. The bisector of a line is de- 



fined as the middle points of the line profil e taking points 
of equal intensity at both sides o f the line l|Toner fc Gravl 
Il98a iLopez-SantiagQ eral]l2003|) . 

The Balmer lines are broadened during the early-flare 
phases and the broadening decreases along the flare evolu- 
tion (see Fig.[7|). Although the emission rises in both wings 
of the Balmer lines, the observed enhancement and broad- 
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Table 7. Width at the base of the chromosphcric lines in 
the quiescent state and at the maximum of flare 2. 



4850 4660 4670 

Wavelength (A) 



3920 



3930 3940 
Wavelength (A) 



Fig. 7. Evolution of the H/3 (left) and Can K (right) hne 
profiles (solid line) during the strongest observed flare 
(flare 2), compared with the quiescent state (dotted line). 



ening are bigger in the red wing. Therefore the Balmer 
lines show a red asymmetry during the detected flares, 
and the largest asymmetry is observed at flare maximum 
(see also the behaviour of the bisectors given in Fig. ISJ . 
Besides, the stronger the flare (flares 2 and 6) the greater 
this asymmetry is. The red asymmetry is also stronger 
in higher members of the Balmer series. We have used a 
two-Gaussian (narrow=N and broad=B) fit to model the 
Balmer lines at flare maxima. The B component appears 
red-shifted with respect to N. For example, at the max- 
imum of flare 2 the distance between both components 
increases from 0.4 A in the case of H/3 up to 1.5 A in Hg. 



jPovle et alJ l)198j^ observed a similar effect during a flare 
on YZ CMi, finding that H7 and US showed symmetrically 
broadened profiles while Hg and Hg showed predominantly 
red-shifted material. They suggested that within an expo- 



Line 


Width at the base of the line 


Quiescent state 


Maximum of flare 2 


H/3 


5±1 


13±1 


H7 


4±1 


6±1 


H5 


4±1 


9±1 


He I 4026 A 


1±1 


1±1 


Can H + He 


4±1 


6±1 


Can K 


3±1 


5±1 


Hg 


3±1 


6±1 


H9 


3±1 


5±1 


Hio 


3±1 


3±1 


Hii 


3±1 


5±1 



sure time as those used in this work (see Table ^) several 
downflows, corresponding to different flare kernels which 
brighten successively one after another, may be present in 
the stellar atmosphere. Every downflow would produce a 
red-shifted contribution to the Balmer lines, which would 
be larger for higher members of the Balmer series (see 
S I3.2|I . In addition, a smaller red asymmetry in the Balmer 
lines is also observed during the quiescent state of AD Leo. 

The Call H & K lines are not usually broadened and 
do not show any asymmetry, excepting the case of the 
strongest flare, where little variations are observed (see 
Fig. Hand HI. 

There are several interpretations for the asymmetries. 
Broad Balmer emission wings have been observed during 
flares produced by very different kinds of stars, showing 
blue or red asymmet rie s llDqvle e t al."l988^ ' Phillips et alJ 
Tflsi lEasonet"^ 119921 iHoudebine et abl Um?i_ 



1988t lEasonet aLl Il992t iHoudebine et a.Ll 1199,^ 
Gunn et alJll994aHAbdul-Aziz et alJll99,'THAbranin et alJ 



T99S iBerdvugina et al f 

iMontes fc RamsevI Il999t 

or II P obvious line asymmetries iHawleT^TPettersenl 
I1991I) . Broadened profiles and red asymmetries ar e 
important constraints on flare models l(Pallavicinilll99(l) . 
They can be attributed to plasma turbulence or mass 
motions in the flare region (see iMontes et alJ Il999t 
iFuhrmeister et al]l2005l and references therein). In solar 
flares, most frequently, a red asymmetry is observed in 
chromosphcric lines. Taking into account that the flares 
reported in this work are like the most typical solar 
flares - non white-light flares - it was expected to find 
red asymmetries like those observed. However, evidences 
of blue asy mmetries hav e been also reported during 
solar flares l|Heinzel et al.l ll994'). Red asymmetries have 
been often interpreted as the result of chromosphcric 
downward condensations (CDC) (see ICanfield et alJ 
I1990I and referen ces therein ) . In f act, recent line pro- 
file calculations llCan et al.l Il993l iHeinzel et all Il994t 
iDing Fan^ll997i) show that a CDC can explain both 
blue and red asymmetries. On the other hand, evidences 
of mass motions have been also reported during stellar 
flares. In particular, a large enhancement in the far blue 
wings of Balmer lines during the impulsive phase of a 



1^ 



iMontes etal J lie 
Lopez- Santiago et al.l 12005 
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Fig. 8. Evolution of the bisector of the H_g (left) and. Call K (right) lines during the strongest flare (flare 2). Note 
that only the minimum value of each bisector is placed on the JD in which its respective spectrum was taken. 



stellar flare was interpre ted as high velocity mass ejection 
ijHoudebine et al.l Il990l) or high velocity chromospheric 
evaporation llCxunn et all Il994bl) . wher eas red asymme- 
tries i n Balmer lines were reported by iHoudebine et all 
^|199.1^ as evidence of CDC. 

Finally, the small red asymmetry observed in the 
Balmer series during the quiescent state can be therefore 
interpreted as multiple CDC happening in the stellar at- 
mosphere. These CDC may be produced by unresolved 
continuous low energy flaring. 

4. Balmer decrement line modeling 

Balmer decrements (flux ratio of higher members to H7) 
have been frequently used to derive plasma densi ties and 
temperatures in the chromosphere of fl are stars llKunke 



IT97ft lCershbei^ri974l lKatsovairi99(i. I.Tevremovic et " 
l|l99j^ developed a procedure (BDFP hereafter) to flt the 
Balmer decrements in order to det ermine some physical 
paran i eters in the flaring plasma. lOarcia- Alvarez et al.l 
l)2002l l2006j) used the BDFP to obtain a detailed trace 
of physical parameters during several flares. 

The BDFP is based on the solution of the radiative 
transfe r equation. It uses the escape pr obability tech- 
nique l)Drakelll98(i iDrake &: Ulrichlll98(]|) and a simpli- 
fied picture of the flaring plasma as a slab of hydrogen 
with an underlying thermal source of radiation which 
causes photoionization. This source represents a deeper 
layer in the stellar atmosphere, which is exposed to ad- 
ditional heating during flares. For a detailed description 
of the physical assumptions made by the BDFP, see 
§ 4.1 of ICarcia- Alvarez et alJ l|20n2l) . The BDFP mini- 
mizes the difference between the observed and calculated 
Balmer decreme nts using a multi-directional search algo- 
rithm l)Torczoni ri991. 1993). This allows us to find the 
best possible solution for the Balmer decrements in a 
four dimensional parameter space, where the parameters 
are: electron temperature (To), electron density (tIc), op- 
tical depth in the Lya line (rLya), and temperature of 
the underlying source or background temperature (Tug). 
Although the temperature of the underlying source could 



be a free parameter in our code, we have fixed a lower limit 
at 2500 K for numerical stability reasons. The best solu- 
tion for the Balmer decrements allows us to calculate the 
effective thickness of the slab of hydrogen plasma and the 
total em ission measure per un it volume (see expressions 
given bv lDrake fc Ulrichlll98nl) . With these quantities we 
are also able to determine the surface area of the emitting 
plasma - volume/thickness - taking into account that the 
volume is the ratio of the observed flux in the H/5 line 
to the total emission measure (also in this line) per unit 
volume. 



4.1. Physical parameters of the observed flares 

We apply the BDFP to the flares observed on AD Leo. 
The Balmer decrements for H7, Hi5, Hg, Hg and Hio have 
been calculated at each flare maximum using the flare- 
only fluxes given in Table ]E\ We do not include higher 
Balmer lines due to both the difficulty in assigning their 
local continuum level and the very low SNR in their spec- 
tral region. 

Fig. 1^ shows the observed and fitted Balmer decre- 
ments at the maximum of the detected fiares. The 
strength of an optical fiare is related to the slope of the 
fit so lution for the Balmer decrements ( Garcia- Alvared 
shallower slopes means stronger flares. The re- 
sults after applying the BDFP are listed in Table |H1 The 
flares 3, 6 and 7 have a relatively high electron den- 
sity (rie > 1 X IQ-'^'' cm~'^) compared with that found for 



the other flares (6 x 10 cm < n. 



< 9 X 10^3 



The flare 3 also has the lowest background tempera- 
ture (Tus ~ 8000 K) which means that the ionization 
balance in this flare is even less radiatively dominated 
than in the other ones (9500 K < < 13500 K). The 
electron temperature of the observed flares ranges from 
12000 K to 24000 K. Taking into account that Hio seems 
to be out of the trend followed by the other Balmer decre- 
ments on flares 2 and 3 (see Fig. O we have also run the 
BDFP rejecting this point, obtaining no significant varia- 
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Table 8. Physical parameters obtained for the observed flares. The area and stellar surface percentage covered by these flares 
are also shown. 



Flare 


Date 


UT 


log TLya 


log Tic 


logTe 


log Tus 


Area 
(x 10^'' cm^) 


Surface 
(%) 


1 


03/04/01 


00:30 


3.71 


13.75 


4.31 


4.13 


1.17 


0.40 


2 


03/04/01 


21:32 


4.62 


13.79 


4.07 


3.99 


6.79 


2.30 


3 


04/04/01 


00:55 


4.18 


14.22 


4.09 


3.90 


1.45 


0.49 


4 


04/04/01 


01:36 


3.78 


13.95 


4.30 


3.98 


0.89 


0.30 


5 


04/04/01 


21:11 


4.54 


13.79 


4.12 


4.02 


0.76 


0.26 


6 


04/04/01 


22:08 


4.43 


14.12 


4.07 


4.04 


3.14 


1.07 


7 


04/04/01 


23:41 


4.41 


14.39 


4.09 


4.02 


2.13 


0.72 


8 


05/04/01 


01:18 


4.15 


13.78 


4.17 


3.98 


2.12 


0.72 


9 


05/04/01 


02:25 


3.60 


13.88 


4.38 


4.09 


0.36 


0.12 



tions on the results (the logarithmic values of the physical 
parameters given in Table [S] change a mean of 5 %). 

In addition, we have used the fit solutions for cal- 
culating the surface area covered by the flaring plasma. 
The values, which are between 3.6 x 10^^ cm^ and 
6.8 X 10^^ cm^, are shown in Table |H1 We find that all 
the observed flares cover less than 2.3 % of the projected 
stellar surface (ttR^). 

All the obtained physical parameters are consistent 
with previously derived values for stellar flares. These 
flares are also similar in size and strength to those 
analyzed in previous works based on the same eg de 
l|Garcia-Alvarezll2003t iGarcia- Alvarez et aljEool |2006|) . 

5. Discussion and conclusions 

The star AD Leo has been monitored using high tem- 
poral resolution during 4 nights (2-5 April 2001). 
More than 600 intermediate resolution spectra have 
been analyzed in the optical wavelength range. 
Although large variations have not been observed, 
we have found frequent short (duration between 
14 ± 1 and 31 ± 3 min) and weak (released energy in 
H/3 between 1.5 x 10^^ and 8.2 x lO^^ erg) flares. Most 
of these fla res are even weaker than those observed on 
AD Leo bv iHawlev et all l)2003l) . All the detected flares 
have been inferred from the variation in the EW of 
different chromospheric lines, which has been measured 
using a very accurate method. Given that the continuum 
does not change during these events, they can be classified 
as non white-light flares, which so far represent the kind 
of flares least known on stars. 

The observed flare activity is > 0.71 flares/hour, 
which is slightly larger than the values that other authors 
obtained for this star by using photometry. As far as we 
know, it is the first time that such a high flare frequency 
is inferred from the variation of the chromospheric lines. 
In fact, AD Leo seems to be continuously flaring since 
we have detected 14 moderate flares and a great num- 
ber of smaller events that appear superimposed on them. 
Given that the energy distribution of flares was found to 



be a power law ijDatlowe et al .ll974tlLin" et al .119841) of the 

form dN/dE = kE~°' - where dN is the number of flares 
(per unit time) with a total energy (thermal or radiated) 
in the interval [E, E + dE], and a is greater than - 
we can expect very weak flares occurring even more fre- 
quently than the observed ones. Therefore our results can 
be interpreted as an additional evidence of the important 
role that flares can play as heating ajgents of the outer at- 
mospheric stellar lay ers (Giidcl 1997: Audard et al. 2000; 
Ka,shva,T^ et alJEfini iGiidel et al. .2003; Arzner Giideli 
20041) . 

A total of 14 flares have been studied in detail. The 
Balmer lines allow to distinguish two different morpholo- 
gies: some flares show a gradual decay much longer than 
the impulsive phase while another ones are less impulsive. 
This resembles the two main ty pes of solar flares (erup- 
tive and confined) described bv IPallavicini et "all l)l977l) . 
However, the Can H & K lines always show the same be- 
haviour and their evolution is less impulsive than that 
found for the Balmer lines. The two maxima and/or weak 
peaks, observed sometimes during the detected flares, 
can be interpreted as the succession of different mag- 
netic reconnection processes, which would be originated 
as consequen ce of the disturbance pro duced by the orig- 
inal flare, as IKopp fc Pneuina 3 lll976l) suggested for the 
Sun. The Can H & K lines reach the flare maximum after 
the Balmer series. It seems that the time delay is higher 
when the impulsive phase of the Balmer lines is shorter. 
The moment in which a line reaches its maximum is re- 
lated to the temperature that characterizes the formation 
of the line and, therefore, is also related to the height 
where the line is formed. During the detected flares, the 
relative increase observed in the emission of the Balmer 
lines is greater for lower wavelengths. We have also found 
that, not only the detected flares are like the most typical 
solar flares (non white- light flares), the asymmetries ob- 
served in chromospheric lines (red asymmetries) are also 
like the most frequent asymmetries observed in the Sun 
during flares. The detected broad Balmer emission wings 
and red asymmetries can be attributed to plasma turbu- 
lence, mass motions or CDC. The Can H & K lines seem 
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Fig. 9. Observed Balmer decrements (squares) and optimum computed fits (solid line) at the maximum of the observed flares. 
The pseudo-quiescent state for each fiare, as defined in is plotted as reference (dotted line). Note that the Balmer decrement 
for Hio is out of the plotted region in fiare 3, but it has been taken into account for doing the fit. 



to be less affected by flare events and their broadening is 
negligible. A small red asymmetry in the Balmer series is 
also observed during the quiescent state, which could be 
interpreted as multiple CDC probably due to continuous 
flaring of very low energy. 

In addition, we have used the Balmer decrements as 
a tracer of physical parameters during f lares by using the 
model developed bv Ijevremovic et al] l)l998|) . The phys- 
ical parameters of the flaring plasma (electron density. 



electron temperature, optical thickness and temperature 
of the underlying source), as well as the covered stellar 
surface, have been obtained. The electron densities found 
for the analyzed flares (6 x 10"^^ cm~'^ - 2 x 10^'* cm~^) 
are in general agreement with those that other authors 
find by u sing semi-empirical chromosphcric modeling 
l)Hawlev fc Fisher 1992a b; Mauas & Falchi 1996). The 
electron temperatures are between 12000 K and 24000 K. 
The temperature of the background source ranges from 
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Fig. 10. Flare duration, stellar surface and temperature 
of the underlying source at flare maximum vs. energy re- 
leased in H/3 during the detected flares. 



8000 K to 13500 K. All the obtained physical parameters 
are consistent with previously derived values for stellar 
flares. The areas - no larger than 2.3 % of the projected 
stellar surface - are c omparable with the size of other so- 
lar and stellar flares llTandberg-Hanssen fc EmslielllflSSt 
lGarcia-Alvare3l2003t ICarcfa- Alvarez et al.ll2002l |2006|) . 

We have also analyzed the relationships between the 
flare parameters. The released energy is correlated with 
the flare duration and the area covered by the flaring 
plasma, but not with the temperature of the underlying 
source (see Fig. . The se results are in gen eral agree- 
ment with those found bv lHawlev et alJ l|2003|) using only 
4 flares and a different method for obtaining the physi- 
cal parameters. We have found a clear relation between 
the released energy and the flux at flare maximum (see 
Fig. 111(1 . Besides, the higher the flux, the longer the flare. 
This flux is also correlated with the area covered by the 
emitting plasma at flare maximum, but not with the tem- 
perature of the underlying source fFig. Ill|) . No correla- 
tions between the area, temperature and duration have 
been found: i.e. these parameters seem to be independent. 
The magnetic geometry seems to be a very important fac- 
tor in flares: firstly, the flare duration can be related t o 
the loop length, as in X-rays (see review bv iReal jl2002l) . 
in the sense that a fast decay implies a short loop and a 
slow decay a large loop; and secondly, the flare area can 
be related to the loop width, because the surface covered 
by the Balmer emitting plasma would be the projected 
area of the flaring loops. On the other hand, the tem- 
perature of the underlying source could be related to the 
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Fig. 11. As Fig. ^1 but vs. the flux emitted in H/3 at flare 
maximum. The top panel has been added to show the 
relation between this flux and the energy released in H/3. 



depth of the layer reached by the flare accelerated parti- 
cles, which is also related to the energy released through 
magnetic reconnection. For explaining the fact that the 
temperature of the underlying source seems to be well- 
defined and independent of the fiare energy (see Fig. [TUl 
and lll|) . we suggest that larger energies may imply more 
energetic particles and therefore larger depths, but not 
necessarily higher temperatures. 

This work has extended the current number of stellar 
flares analyzed using high temporal resolution and good 
quality spectroscopic observations. However, there is still 
need for new data of this type to trace the behaviour of 
the physical parameters throughout these events. This will 
help us to understand the nature of fiares on dMe stars. 
Finally, higher spectral resolution would be also desirable 
in order to study the changes that take place in stellar 
atmospheres during flares. 
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